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Abstract 

The movement of a Brownian particle in a fluctuating two-state periodic potential is investigated. 
At high transition rate, we use a perturbation method to obtain the analytical solution of the 
model. It is found that the net current is a peaked function of thermal noise, barrier height and 
the fluctuation ratio between the two states. The thermal noise may facilitate the directed motion 
at a finite intensity. The asymmetry parameter of the potential is sensitive to the direction of the 
net current. 
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I. INTRODUCTION 



Much of the interest in non-equilibrium-induced transport processes has been on the 



stochastically driven ratchets 



>n-equiii brium-mcu. 



The noise-induced ratchet has recently 
attracted considerable attentions. It is related to the symmetry breaking and comes from the 
desire for an explaination of directional transport in biological systems. Several models have 

and the asymmetric 



been, for example, proposed to describe the muscle contraction [3, Q] 
polymerization of action filaments responsible for cell mobility 

The focus of research has been on the noise-induced unidirectional motion over the last 
decade. In these systems directed-Brownian-motion of particles is generated by nonequilib- 
rium noise in the absence of any net macroscopic forces and potential gradients. Ratchets 
have been proposed to model the unidirectional motion due to the zero-mean nonequilib 

n n " 

rium fluctuation. Typical examples are rocking ratchets |10l . I 111 . 
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1 131] .flashing ratchets 

diffusion ratchets M and correlation ratchets M- Ghosh et. al.Q Q devefoped 
some analytical solutions of the current variation with the noise in the ratchet. In all these 
studies the potential is taken to be asymmetric in space. It has also been shown that a uni- 
directional current can appear for spatially symmetric potentials. For the case of spatially 
symmetric potential, an external random force should be either temporally asymmetric or 
spatially-dependent. 

The previous works are limited to case of single potential. The present work extends the 
analysis to the case of two potentials in flashing thermal ratchet: one constant potential 
and one periodic in space and constant in time. No external driving forces are required to 
induce unidirectional current in these ratchets of two potentials. The emphasis is on the 
current as the function of noise and other system parameters. This is achieved by using a 
perturbation method to solve two coupled Simoluchowsky equations. 



II. FLASHING MOTOR 



Consider the flashing motor (flashing ratchet) model aiming for describing the spatially 
unidirectional motion along x-direction in Fig. 1 of a Brownian particle due to the two 
potentials (states): one constant and the other spatially-periodic. This model was initially 
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proposed in an attempt of describing molecular motor in biological systems [191 ]. 

The rate of fluctuation between the two potential states is governed by two rate constants 
k\ and k 2 , respectively. Here the former is the rate from State 1 to State 2 and the latter is 
the rate from State 2 to State 1. While the particle diffuses freely at State 2, it is localized 
near a local minimum at State 1. The particle motion satisfies the dimensionless equation 
of motion 

where is the Brownian random force, x the position of the particle, t the time, Vi(x) 

the potential, the subscript i stands for the state which can take value of 1 or 2. The prob- 
ability densities Pi(x, = 1,2) of state 1 and 2 are govern by two coupled Simoluchowsky 
equations 

dP2(x,t) „ &P 2 {x,t) dJ 2 

Qj. = D —fa2 + hPi(x,t) - k 2 P 2 (x,t) = - — , (3) 

where f(x) = —V 1 (x), the prime stands for the derivative with respect to x, Ji, J 2 are 
probability densities of current, D the diffusivity, k% Boltzmann constant, T the absolute 
temperature. Here x, t, k±, k 2 , D, ksT are all dimensionless. Also, 

. ^(x — mL), mL < x < mL + A; 

V 1 (x) = { xK h ~ (4) 

■^[-x + (m + 1)L], mL + \<x<(m + 1)L, 

where m = 0, 1, 2, .... At steady state such that ^ = and ^ = 0, Eqs (2) and (3) lead 



to the net current 



where 



and 



J=-D^&+D 1 f{x)P 1 {z), (5) 
D[P"{x) - P'i{x)\ + k[(l + /i)Pi(x) - tiP(x)] = 0, (6) 



J=J, + J 2 , P(x) = P 1 (x) + P 2 (x), D 1 = -f-, (7) 

k B T 



k = ki, \i — k 2 /k. (8) 
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FIG. 1: Two-state model with regular ratchets: V\{x) is spatially-periodic sawtooth of period L 
and barrier height Vo- A is an asymmetry parameter (Vi(x)is symmetric when A = L/2); V2(x) is 
a constant potential. 

III. ANALYTICAL SOLUTION 

When the fluctuation is at high rate such that k » 1, we can expand P{x), P\{x) and 
J in power series of a small parameters fc -1 j^. 

oo oo oo 

P{x) = ]T k- n Pn (x), P x {x) = £ k- n Pln (x), J = £ k~ n Jn . (9) 

71=0 n=0 n=0 

The coefficients p n , p\ n and j n can be obtained by substituting Eq.(9) into Eqs (5) and 
(6) and equating coefficients of k~ n , 

Po(x) - ^^ D /(^)Po(^) = (10) 



Pio{x) = — 7— PoO), (11) 

1 + /U 



Dp^a) + ^-f(x) Pn (x) = 3n + G n ^(x),n = 1, 2, 3, (12) 
G n(a;) = Y^-f{x)\p" n {x) -pl n {x)],n = 0,1,2,.... (13) 

-L ~T~ \-k 
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Under the periodicity conditions 



Pn(x + L) =p n (x), 71 = 0,1, 2,... (14) 

and the normalization of the distribution p(x) over the period L, 

[ p n (x)dx = 5 0n ,n = 0, 1,2, (15) 
Jo 

we can obtain all coefficients of p n , p\ n and j n . Since our attention is mainly on the current 
J, we only list j n (x) here 

Jo = 

_ £ G n - X {x)U-\x)dx 
3n ~ ! L U~\x)dx ' n " i ' 2 '- W 



where 



In particular, 



U(x)=e^[-jJ^ 5 V 1 (x)]. (17) 
Jl (1 + J Q L U(x)dxJ L U- 1 (x)dx { ' 



Therefore, to the first-order approximation, 



j ~ o + = ^1 Jof 3 M d * (19) 

k{l + fi) 4 D j L U{x)dxj L U- 1 {x)dx 



After substituting f(x) and U(x), we have 

^ 4 £V 5 (2A - L) 
(1 + fi) 6 f3 5 kLX 2 {L - \) 2 {eT^ Vo + e"(TT7iF ,/0 - 2) 

where (3 = ksT. By letting |^ = with 7 = ^, we have 



<*» " T^) + H- j^) exp[^] - 10, = 0, (21) 



which leads to the optimum 7 for the maximum J (J moXj7 ). 



By letting |^ = 0, similarly, we have the optimum \x for the maximum J (J ma x,fi), 



- 2/ i 2 + (2--)/i + 4)]exp[ » ] + [ _2 At 2 + (2 + l) / i + 4)]exp[ / ]+4(/x 2 -//-2) = 0. 

(22) 
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FIG. 2: Plot of the optimum 7 vs the optimum fi for the maximum J [Vq = 5.0, D = 1.0, L = 2.0, 
fc = 100.0, A = 0.5). 




FIG. 3: Plot of the maximum J shown in Eqns. (21) and (22). (a) Plot of the maximum J with 
the optimum 7 vs [i. (b) Plot of the maximum J with the optimum fx vs 7 (Vb = 5.0, D = 1.0, 
L = 2.0, jb = 100.0, A = 0.5). 



IV. RESULTS AND DISCUSSION 

Equation (20) indicates that the direction of the net current is determined by the asym- 
metry parameter A. When < A < 1, the current is positive, the current is negative when 
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1 < A < 2. There is no current at A = 1. 

Figure 2 shows the solution of the Eqns (21) and (22). The dot line gives the optimum 
7 for the maximum J (Eq. (21)) and the solid line gives the optimum \x for the maximum 
J (Eq. (22)). It is easy to find that the dot line meets the solid line at the point A 
((J>opt,Jopt), which indicates that one can obtain the maximum J for both the optimum 7 and 
the optimum /x at the same time. The corresponding J maX)7 vs /1 and J ma x,^ vs 7 are shown 
in Fig. 3a and Fig. 3b, respectively. From Fig. 3a, we can find J ma x,-y as the function of /1 
have a maximum value, at which the \x and 7 are optimal, namely, the solid line will meet 
the dot line as shown in Fig. 2. The similar results can also be obtained in Fig. 3b. 
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FIG. 4: Dimensionless probability current J vs thermal noise strength k B T for different values of 
asymmetric parameters (Vq = 5.0, D = 1.0, L = 2.0, k = 100.0, \i = 1.0). 

Figure 4 shows the variation of the net current J with the thermal noise intensity ksT. 
The curve is observed to be bell-shaped, a feature of resonance. The current reversal appears 
at A = 1 at which the potential Vi(x) is symmetry. When kBT — > 0, J tends to zero for 
all values of A. Therefore, there are no transitions out of the wells when the thermal noise 
vanishes. When k B T — > 00 so that the thermal noise is very large, the ratchet effect also 
disappear. The current \J\ has a maximum value for fixed value A at certain value of ksT. 
By Eq. (21), the optimized value is k B T = 0.5073(7 — 0.1015, Vq = 5.0). The maximum 
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Jmax = 0.0949 at A = 0.5. Therefore, certain thermal noise can induce a large current | J 
while the thermal noise blocks the unidirectional motion in general. 
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FIG. 5: Dimensionless probability current J vs barrier height Vq for different values of A (fceT = 0.5, 
D = 1.0, L = 2.0, k = 100.0, n = 1.0). 

Figure 5 shows the net current as the function of barrier height V . We again observe the 
current reversals at A = 1.0. When the barrier height Vq is small. The effect of the ratchet 
is also small; the thermal noise effect is domaint so that the net current disappear. When 
the barrier height V is large, on the other hand, the particle can not pass the barrier. It 
can only diffuse at State 1 so that the net current is also very small. Therefore, there is 
an optimized value of V (4.9281) at which J takes its maximum v&\ue(J max = 0.0949), for 
example, A = 0.5. 

Figure 6 shows the current as the function of fi. When fi — > 0, /C2 tends to zero. The 
attraction by State 2 is too small such that the particle is always staying at State 1. The 
ratchet reduces to one-state ratchet without external force. Therefore, no current exits. 
When /i » 1, the attraction by State 1 becomes too small so that the particle can only 
diffuse at State 2. The net current becomes to zero again. Hence, there exits an optimized 
value of fi (0.3991) at which J takes its maximum value(J max = 0.1954) as shown in Eq. 
(22). 
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FIG. 6: Dimensionless probability current J vs the ratio of the two state transition rate fi for 
different values of A {k B T = 0.5, V = 5.0, D = 1.0, L = 2.0, k = 100.0). 

V. CONCLUDING REMARKS 

Two coupled Simoluchowsky equations are solved by a perturbation method to obtain 
the net current. The current is peaked function of thermal noise k B T, barrier height Vq and 
the fluctuation ratio /x between the two states. It is positive for < A < 1 and negative for 
1 < A < 2. Therefore, the current reverses its direction at A = 1.0 (symmetric potential). 
When the thermal noise is small, the particle can not pass the barrier such that the current 
J tends to zero. When the thermal noise is too large, the ratchet effect disappear so that 
J tends to zero, also. There is an optimized value of thermal noise at which J takes its 
maximum value. For the case of Vq — > 0, the thermal noise is dormant and the current 
disappears. When Vo — > oo, on the other hand, the particle can not pass the barrier. When 
/j, — > 0, the attraction from State 2 is too small and the particle is always at State 1. The 
ratchet reduces to one-state ratchet without external force. Therefore, no currents occur. 
When /i is very large, similarly, the attraction from State 1 is too small, the particle can 
only stay at State 2, and the current also tends to zero. There exits optimized values of 
ksT,Vo and /i at which the current takes its maximum value. 

Here, the thermal noise can facilitate the directed motion of the Brownian particles. 
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This differs from the one prevalent in the literature that the thermal noise always de- 
stroyed the directed motion. The noise-induced transport is associated with the breaking 
of either reflection symmetry of spatially periodic system or statistical symmetry of 
temporal nonequilibrium fluctuations characterized by multitime correlation functions. 
The symmetry-breaking driven transport can be used to explain the directed motion of 
macromolecules in biological cell and to construct well-controlled devices of high resolution 
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for separation of macro-particles and micro-particles 
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